Berry Phase and Pseudospin Winding Number in Bilayer Graphene 
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Ever since the novel quantum Hall effect in bilayer graphene was discovered, and explained by a 
Berry phase of 2n [K. S. Novoselov et at, "Unconventional quantum Hall effect and Berry's phase 
of 2n in bilayer graphene", Nature Phys. 2, 177 (2006)], it has been widely accepted that the 
low-energy electronic wavefunction in this system is described by a non-trivial Berry phase of 2n, 
different from the zero phase of a conventional two-dimensional electron gas. Here, we show that 
(i) the relevant Berry phase for bilayer graphene is not different from that for a conventional two- 
dimensional electron gas (as expected, given that Berry phase is only meaningful modulo 2n) and 
that (ii) what is actually observed in the quantum Hall measurements is not the absolute value of 
the Berry phase but the pseudospin winding number. 
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I. INTRODUCTION 

A conventional two-dimensional electron gas (2DEG) 
under a high magnetic field and low temperature shows a 
quantized Hall conductance a xy — n e 2 /h per spin degree 
of freedom where n is an integer, e the charge of an elec- 
tron, and h Planck's constant [1[. For graphene, it has 
been predicted that the Hall conductance per spin and 
valley degrees of freedom is a xy = (n+l/2)e 2 /h 03, 
and this prediction has then been proven by experiments 
performed on mechanically exfoliated samples [B|, 
Such half-integer quantum Hall effect is the most decisive 
evidence that the sample is actually a single atomic layer 
of carbon atoms effectively decoupled from the substrate, 
and is a direct manifestation of a non-trivial Berry phase 
of 7r in the electron wavefunction d@]. (This half-inte gcr 
quantum Hall effect is different from the fractional quan- 
tum Hall effect in graphene 

[Mo 

Half-integer quan- 
tum Hall effect has been observed in graphene epitaxially 
grown on the carbon-rich surface [111 ] and the silicon-rich 
surface 12h14|] of silicon carbide or in graphene grown by 



chemical vapor deposition |15j . 




FIG. 1: (a) Schematic representation of graphene and bilayer 
graphene. For graphene, the carbon atoms belong only to the 
two A and B sublattices. (b) The Brillouin zone of graphene 
and bilayer graphene. The positions of the K and K' points 
are (|f-,0) and (— §J,0), respectively, where a is the lattice 
parameter. 



When a second layer is added — thus forming bilayer 
graphene [Fig. 1(a)] — the electronic structure changes 
dramatically. The dispersion for low-energy quasipar- 
ticles in graphene is linear [Fig. 2(a)], whereas in bi- 
layer graphene it becomes quadratic [Fig. 2(b)]. The 
quantum Hall conductance in bilayer graphene has been 
theoretically predicted [l6| . and experimentally verified 
(again by using mechanically exfoliated samples) [13] to 
be <j xy = n e 2 jh per spin and valley degrees of freedom. 
This is very similar to the case of a conventional 2DEG, 
but, very significantly, with no plateau at n = 0. This 
novel quantum Hall effect in bilayer graphene has been 
explained by the appearance of a non-trivial Berry phase 
of 2-7T in the electronic wavefunction 16, I17JI, different 



from the case of a conventional 2DEG. These pioneering 
studies 16J, [17J have significantly influenced the commu- 
nity working on the novel physics of graphene nanostruc- 
tures, and the concept of a non-trivial Berry phase in 
the electronic wavefunction of bilayer graphene has been 
commonly accepted thereafter. 

In this paper, we show that the Berry phase in the 
electronic wavefunction of bilayer graphene is equivalent 
to that of a conventional 2DEG (i. e. , Berry phase of 27r, 
47r, etc., are not non-trivial but are all equivalent to 0) 
by adopting simple arguments on the gauge freedom of 
electronic wavefunctions. We then reinterpret the origi- 
nal experimental results |l7J : although the quantum Hall 
effect in bilayer graphene is qualitatively different from 
that in a conventional 2DEG, what we need for an ex- 
planation of that novel phenomenon is not a non-trivial 
Berry phase of 2-7r but the concept of pseudospin winding 
number, i. e. , how many times the pseudospin vector - 
the vector representing the relative phase on the two dif- 
ferent sublattices of carbon atoms [18| [see lower panels 
of Figs. 2(a) and 2(b), and Eqs. Q and (JE])] - rotates 
when the electronic wavevector undergoes one full rota- 
tion around the Dirac point [Figs. 1(b) and 3]. We note 
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FIG. 2: (a) Upper panel: Energy dispersions for low-energy electrons in graphene. Lower panel: Pseudospin distribution for 
electronic eigenstates in graphene on an equi-energy contour specified by the dashed curve in the upper panel. The arrows 
represent the direction of pseudospin. Here, we consider the electronic states with wavevectors near the Dirac point K. (b) 
and (c): Similar quantities as in (a) for bilayer graphene and a conventional 2DEG, respectively. Note that the electrons in a 
2DEG have an energy minimum at the center of the Brillouin zone and do not have a pseudospin degree of freedom. 
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FIG. 3: Schematics of electronic eigenstates in graphene 
or multi-layer graphene on an equi-energy contour with a 
wavevector around the Dirac point K. We define q = k K 
where k is the 2D Bloch wavevector. If we start from the state 
|0), we end up with the state \N) after a complete counter- 
clockwise rotation along the contour (|JV) = |0)). Here, the 
evolving parameter is q. 



that this pseudospin winding number is exactly the same 
quantity as the 'degree of chirality' J of previous stud- 



ies 



III 17 1 . On the other hand, the term 'chirality' - in 
connection with neutrino physics - has been used to refer 
to the inner product between the two-dimensional (2D) 
Bloch wavevector from the Dirac point and the pseu- 
dospin vector [lij]; hence, for graphene, the chiralities for 
electronic states with wavevector close to the K point in 
the upper band and in the lower band would be +1 and 
— 1, respectively, whereas the degree of chirality is +1 
regardless of the band index. To avoid this confusion, 
we suggest the term 'pseudospin winding number' here. 
Also, pseudospin winding number has been defined in the 
form of a more general topological invariant 



II. BERRY PHASE 

For electrons in a periodic system, the Berry phase is 
a phase acquired by the wavefunction over the course of 
a cyclic evolution of the Hamiltonian - such concept has 
had a fundamental role in developing the modern theory 
of polarization 23, 24 1 and of magnetization [2ol |26|. 
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The evolving parameter that we consider here is the 
2D Bloch wavevector k = (k x ,k y ). The Berry phase T 
described in the wavevector evolution of Fig. 3 is given 
by 



N-l 

T = -i Jim V log + 1) 

3=0 



where 



\N) = |0) . 



(1) 



(2) 



Here, k is determined by the state index j. 

We first show that Berry phase [defined in Eq. ([I])] is 
in general meaningful only modulo 27r. After multiply- 
ing each of the wavefunctions by a position-independent 
overall phase factor e l 6j , we obtain another set of wave- 
functions which satisfies the equation of motion 



\j)' = e* e ' \j) 



(3) 



In order to have a well-defined set of wavefunctions for 
Berry phase evaluation [Eq. (|5J)], we have to impose 



On — do = 2nm 



(4) 



where m is an arbitrary integer. Now using Eqs. (HJ) , d3J) , 
and (|4|) , the Berry phase T' for the new set of wavefunc- 
tions [Eq. ©] is 



N-l 



V = -i lim ^logO'lj + l}' 



3=0 
N-l 



r+ lim 0,0 



j=0 



T + 27rm, 



(5) 



i. e. , Berry phase can be defined modulo 2ir only. In the 
following, we paraphrase this general argument using the 
language of graphene. 

For graphene, let us consider the case where k is very 
close to the Dirac point K [Fig. 1(b)], and define q = 
k — K (|q| <C |K|). Then, if we use a basis set composed 
of the Bloch sums of p z orbitals localized on the two 
sublattices A and B [Fig. 1 (a)] , the effective Hamiltonian 
reads El 



ffmono(q) =flV q 



exp(— i9q) 
exp(i# q ) 



(G) 



where vq is the band velocity and 6* q the angle between 
q and the +k x direction. The energy eigenvalue and 
wavefunction of Eq. (J6]) are given by E™° no = hvo s |q| 
and 



(7) 



respectively, where s = ±1 is the band index. As shown 



before [27] . the Berry phase r mono of the electron wave- 



function in graphene when the Bloch wavevector under- 
goes a full rotation around the Dirac point (Fig. 3) can 
be obtained from Eq. (fTJ): 



N-l 

—i lim } log 

j=0 

[-21T 

—I 



1 + cxp [i(9 j+ x - Qj)] 



d6- 
2 



7T: 



(8) 



where we have substituted q by the corresponding state 
index j (j = 0, 1, 2, • • • , N — 1) as shown in Fig. 3. 

Now we show that Eq. ^ is not affected by an arbi- 
trary phase factor of each electronic wavefunction. Sup- 
pose that we use a new form for an electron wavefunction 
in graphene 



)' - 



71 v sel9]/2 



(9) 



for j = 0, 1, N - 1, and of course |^™°)' = H'T.o™)' 
(Fig. 3). Then the Berry phase rj„ oao is given by 

rl„ = -i Km {log<^!_ il^o ono )' 

( e -i{e 1+ x-6j)/2 + e i(e j+1 -ej)/2 

+ E l0 ? 

J=0 



2tt 



(10) 

i. e. , we obtain the same result as in Eq. ([8]). The phys- 
ical reason for this invariance of the Berry phase with 
respect to an arbitrary gauge phase of each electronic 
state evaluated from Eq. ((T|) is that each state appears 
twice (once as a bra state and once as a ket state), thus 
cancelling out any arbitrary phase in each wavefunction. 
Since the Berry phase obtained from Eq. (fTJ) is not af- 
fected by any arbitrary phase of each state \j), we can 
choose for convenience a gauge such that the state varies 
smoothly over the course of parameter evolution. 

Similarly, for bilayer graphene, if we use a basis 
set composed of Bloch sums of localized Wannier-likc 
p z orbitals on each of the two sublattices A and B' 
[Fig. 1(a)] [28], and consider only the nearest- neighbor 
intra-layer hopping and vertical interlayer hopping (char- 
acterized by the hopping integral t±) processes, the effec- 
tive Hamiltonian of low-energy electronic states in bilayer 
graphene becomes [17 1 



#bi(q) 



h 2 q 2 
2m* 



exp(-2i6> q ) 
exp(2ifl q ) 



(11) 
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where m* = 3 a 2 i^/8uQ is the effective mass and a the 
lattice parameter. The energy eigenvalue and wavefunc- 
tion of Eq. (TTT|) are given by E^ 1 = sh 2 q 2 /2m* and 



1 



1 



(12) 



respectively (s = ±1). The Berry phase Tbi of the elec- 
tron wavefunction in bilayer graphene can again be ob- 
tained from Eq. ([1]): 



N-l 



r bi 



lim > log 

7_S.~, Z-^ 



Jo 
2tt, 



1 + exp [2i((9 i+ i - 6j)] 



(13) 



which is the result of previous studies 16L 17|. 

Now we consider another form for the wavefunction 



(14) 



one is 



whose only difference from the original 
the position-independent overall phase (IV'sq) = 
e ~i9q Because the two wavefunctions are dif- 

ferent only in the overall coefficient, they both are per- 
fectly good solutions of the effective Hamiltonian of bi- 
layer graphene [Eq. (|lip]]. Moreover, both the original 
[Eq. (|12l) ] and new [Eq. (fl4)) ] wavefunctions satisfy the 
gauge condition that the wavefunction varies smoothly 
in the course of cyclic evolution; hence, in the evaluation 
of the Berry phase, one does not have to consider any 
branch-cut as we did in the derivation of Eq. (ITU1) . The 



Berry phase T' hi for jV'sq)' evaluated from Eq. ([I]) is 



Tbi = 0; 



(15) 



in other words, the Berry phase of the electronic wave- 
function of bilayer graphene is the same as that of a con- 
ventional 2DEG. 

In fact, the Berry phase of an electronic wavefunction 
in bilayer graphene evaluated from Eq. (fTJ) is 2m7T, with 
to an arbitrary integer. Or, in other words, any Berry 
phase of 2 to tt is equivalent to a trivial Berry phase of 0. 
[Similarly, this Berry phase for graphene obtained from 
Eq. (TT]) is (2m + 1) tt, with m an arbitrary integer.] This 
statement is related also to the well known application 
of Berry-phase physics to the modern theory of polariza- 
tion 23|, |24| , according to which the electrical polariza- 



tion in a periodic system can only be determined modulo 
eR/!l, where R is a lattice vector and f2 the unit cell 
volume. 

Our claim that the Berry phase in bilayer graphene 
is not non-trivial and is equivalent to that of a conven- 
tional 2DEG remains valid even if we consider a more 



complex model Hamiltonian. If, for example, second- 
nearest-neighbor inter-layer hopping processes are con- 
sidered in the tight-binding model of bilayer graphene, 
the low-ener gy electronic states can be described by four 
Dirac cones [16| . It was suggested that the central Dirac 
cone contributes —tt to the Berry phase and each of 
the three satellite Dirac cones contribute tt to the Berry 
phase and hence the total "non-trivial" Berry phase is 
— tt + 3 x tt = 2tt [l6[. According to our discussion, this 
view is not correct, again because a Berry phase of — tt 
is equivalent to that of tt, etc. In par ticular, starting 
from the wavefunctions used in Ref. [16[ , one can define 
a new set of wavefunctions with additional four phase fac- 
tors, e.g. , e~ l9 ' 1 , defined around each of the four Dirac 
cones in a way similar to Eq. ([T4")l . Then, the Berry 
phase around the i th Dirac cone (i = 1,2,3,4) obtained 
by evaluating Eq. ([1]) would be (2mi + l)TT for an arbitrary 
integer m^. Thus, the Berry phase for bilayer graphene is 
of the form Tbi = X)i=i(2 TO i + 1)tt = 2mTT, with m an ar- 
bitrary integer, and cannot uniquely be 2tt as suggested 
in previous studies 

sua. 



III. PSEUDOSPIN WINDING NUMBER 

However, these considerations do not mean that there 
is no qualitative difference between the electronic eigen- 
states in bilayer graphene and those in a conventional 
2DEG. Indeedj_a new phenomenon in bilayer graphene 
was observed [13] and explained [l6[: the step size at 
charge neutrality point in the quantum Hall conductance 
a xy in bilayer graphene is 2e 2 /h per spin and valley de- 
grees of freedom and not e 2 /h as in a conventional 2DEG. 
We reinterpret this result as a direct manifestation of the 
pseudospin winding number, and explain this in the fol- 
lowing. 

For an electronic state whose wavefunction can be de- 
fined as a two-component spinor, as in graphene [Eq. ([7])] 
or in bilayer graphene [Eq. (| 1 2[) ] , a 2D pseudospin vector 
can be defined from the relative phase of the two compo- 
18j |. Figures 2(a) and 2(b) show the pseudospin 
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nents 

vector for electronic wavefunctions in graphene and in 
bilayer graphene, respectively. The pseudospin winding 
number n w is then defined as the number of rotations 
that a pseudospin vector undergoes when the electronic 
wavevector rotates fully one time around the Dirac point: 
it is immediate to see from Figs. 2(a) and 2(b) that n w 
is 1 in graphene and 2 in bilayer graphene. 

The low-energy effective Hamiltonian of a graphene 
multilayer, with n layers and ABC stacking, is [29( 



H n (q) (X q n 



exp(— niOq) 



(16) 



exp(n i 9q) 
Eigenvalues and eigenfunctions of Eq. (|16p are given by 



5 



E™ oc s q n and 



cussions, and Intel Corporation for support (CHP). 



(17) 



respectively (s = ±1). Thus, the pseudospin winding 
number n w of the eigenstate in Eq. (|16[) is n. 

It has been shown (Supplementary Information of 
Ref. |l7[) that a system described by Eq. (TTB1) has a quan- 
tum Hall conductance step at the charge neutrality point 
of size n e 2 jh per spin and valley degrees of freedom. Our 
interpretation is that what we learn from a measured step 
size n e 2 /h per spin and valley degrees of freedom |17j is 
not the "absolute" value of the Berry phase (r = nit) 
but the pseudospin winding number {n w = n); We also 
note that the quantum Hall conductance step, which has 
been correctly obtained in Ref. [l7| from Eq. (jTHJ) , does 
not need a Berry phase of V = nir for its explanation. 
In other words, the size of the quantum Hall conduc- 
tance step at the charge neutrality point is connected di- 
rectly to the special form of the Hamiltonian in Eq. (fl~6| , 
or, equivalently, to n which can be interpreted, from its 
wavefunction [Eq. (|17l) ]. as the pseudospin winding num- 
ber and not to the absolute value of the Berry phase. (It 
should be stressed that our discussion of the pseudospin 
winding number is confined to systems whose effective 
Hamiltonian is given by Eq. (1161) .) 

On the other hand, the Berry phase still determines 
the shift of the plateaus, with integer and half-integer 
quantum Hall effects for T = and T — 7r, respec- 
tively 0, [|1 [27} . The former and the latter correspond to 
multi-layer graphene with even and odd number of lay- 
ers, respectively (29j . It has been recently shown within a 
semiclassical theory that the so-called topological part of 
the Berry phase, which is determined by the pseudospin 
winding number, is responsible for this shift [301 ] . 



IV. CONCLUSION 

In conclusion, we have shown that the Berry phase 
of the low-energy quasiparticle wavefunction of bilayer 
graphene is the same as that of a conventional two- 
dimensional electron gas, and that the fundamental dif- 
ference between the two systems is the pseudospin wind- 
ing number. Our findings have a broader implication 
than the cases discussed here. For example, the interpre- 
tation of the recent angle-resolved photoemission spec- 
troscopy on bilayer graphene 3l|, [32( and that of the re- 
cent quantum Hall experiments on trilayer graphene 33 1 
should be revisited in view of the present discussion. 
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